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Thermally activated escape rates of uniaxial spin systems with transverse field: 

Uniaxial crossovers 
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Classical escape rates of uniaxial spin systems are characterized by a prefactor differing from and 
much smaller than that of the particle problem, since the maximum of the spin energy is attained 
everywhere on the line of constant latitude: 9 — const, < ip < 2n. If a transverse field is applied, a 
saddle point of the energy is formed, and high, moderate, and low damping regimes (similar to those 
for particles) appear. Here we present the first analytical and numerical study of crossovers between 
the uniaxial and other regimes for spin systems. It is shown that there is one HD-Uniaxial crossover, 
whereas at low damping the uniaxial and LD regimes are separated by two crossovers. 

PACS numbers: 05.40.-a, 75.50.Tt 



The study of thermal activation escape rates of fine 
magnetic particles, which are usually modelled as clas- 
sical spins with predominantly uniaxial anisotropy, may 
be traced from the early predictions of Neel jl[] through 
the first theoretical treatments of Brown [|]|| to the re- 
cent experiments of Wernsdorfer et al Q on individual 
magnetic particles of controlled form. These experiments 
allow one for the first time to check the Stoner-Wohlfarth 
angular dependence of the switching field j5j and to make 
a comparison [^j?]] with existing theories where the energy 
barrier is reduced by applying a magnetic field. The theo- 
ries checked are those for the intermediate-to-high damp- 
ing (IHD) case as well as for the low-damping 
(LD) case @. 

The IHD and LD limits for spins are similar to those for 
the particle problem, which were established by Kramers 
jllj. The most significant difference is that for spins 
in the HD limit the prefactor To in the escape rate 
r = r exp(— AU/T) behaves as r cx a, where a is the 
damping constant [if the Landau-Lifshitz (LL) equation 
is used], whereas for particles T oc 1/a. A question 
which has not yet been addressed, both theoretically and 
experimentally, and which is the subject of this Letter, 
is how these three damping-governed regimes merge into 
the single uniaxial regime Q if the field is removed? 

Let us consider the Hamiltonian 



the spins f(9, p, t), which follows from the stochastic LL 
equation, reads 
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where x = cos 9, 
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To solve the FPE at low temperatures, T <C AC/, we 
represent / as /(N,i) = / eq (N) 5 (N, t), where / cq (N) = 
Z^ 1 exp[— 7i(N)/T] is the equilibrium distribution func- 
tion. Neglecting the exponentially small g one obtains 
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The function g assumes the values g\ and gi in the wells 
and changes in a narrow region about the top of the bar- 
rier, so that 



g(x, p) = gi + (32 - gi)C(x, <p), 



(•5) 



H = -Kn 2 z - ftnH, 



1, 



where fj,Q = M S V is the magnetic moment and K = KV 
is the uniaxial anisotropy energy of the particle. The 
Fokker-Planck equation for the distribution function of 



(1) where £ assumes the values and 1 in the first and 
second wells, respectively. The numbers of particles in 
the wells satisfy iVi + N2 = 1 and Ni = giNi tCq , where 
A^eq = Zi/Z are the equilibrium values and Z\ and Zi 



are partition functions for each of the wells. 



1 



The change in the particle number's in the 1st well N± 
is due to the flow of particles from 1st to the 2nd well 
through the line x — const (say, the equator x = 0): 
Ni = - J 27r dipj x . Using f cq dH/dtp = -Tdf cq /dip, one 
can integrate the first term of j x by parts. Finally, using 
Eq. (||), one has at the kinetic equations 



Ni = -N 2 = T{N 2 N hcq - N 1 N 2 , e ^, 
the escape rate being 



(6) 



1 { 1 

Zi + Z 2 



7 T 



Mo Jo 



dipe 



—TL/T 



a(l — x ) 



dx dip 



(7) 



where £ satisfies Eq. (jfy with the boundary conditions 
stated in Eq. (||). Also, one can calculate a current 
through a line of constant energy surrounding one of the 
wells, which is appropriate in the LD case. 
For the transverse-field (unbiased) model 



H(x, ip) /T ee e = -o(x 2 + 2/i-v/i - x 2 cos <p), (8) 
where 



K 



<J = — — 



T 



KV 
-jr, 



h = 



2JT 



2JST 



(9) 



Another useful dimensionless variable is £ 
HM S V/T. The function C satisfies 



2crh 



(1 - hcosipjx a- — h Tr - 
\ ox dtp 

2ct \dx* dip 2 ) 



hsmip [dx-- a d^ 



(10) 



where 1 



a; 



1 in the region x <C 1, which is relevant 



in the uniaxial and IHD cases. The result for T can be 
written 



T = — Aexp 

7T 
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Vl-h 2 , (11) 



where AU = H sad - H min , AU /T = tr(l - /i) 2 and 
wi is the ferromagnetic resonance frequency near the 
bottom of the well (the attempt frequency), the factor 
A describes deviations from the transition-state theory 
(TST), and the lack of the factor 2 in the denominator 
is due to the two symmetric wells. 

In the uniaxial case, h — 0, the function £ is inde- 
pendent of ip, and the solution of Eq. ( [To| ) is d(/dx — 
^a/Trexp(-ax 2 ). Eq. (@) then yields @] 



A = 27TO1 



(12) 



Another analytically solvable case is IHD with a pro- 
nounced saddle: T <C Ti. max — Tt sa d, i.e., £ = 2ah ^> 1 
for our model. Here the coefficients in Eq. (|lTJ) for £ can 
be linearized near the saddle point x — ip — shown in 
Fig. [TJ Next one seeks a solution ((x, ip) = Q(u) with 



x (arb. scale) 




FIG. 1. Trajectories of a magnetic particle with uniaxial 
anisotropy and transverse field. Regions 1 and 2 are the po- 
tential wells, O is the saddle point, 3 is the maximum of the 
potential. 

u — x + vip and appropriately chosen v ||[ll) . This form 
of £(x, ip) implies that it changes for a > 1 in a nar- 
row region across the line AOB' in Fig. [j]. Evaluating 
d((u)/du, from Eq. (^) one obtains || 



_ <z(l - 2h) + y/a 2 + 4h(l - h) 
2y/h(l - h) 

in Eq. (|ll|), with the limiting forms 



1, a 2 <^h{l-h) 
\'={ \/i - h(l + a 2 )/a, l-fr<l,a 2 
a/v^, /i < 1, a 2 . 



(13) 



(14) 



Note that for h <C 1 the HD regime is attained already for 
a > y/h <C 1, where the difference between the LL and 
Gilbert equations is still irrelevant. For h < 1/2, Eq. 
([TJ|) monotonically increases with a. For h > 1/2, it has 
a minimum at a = 2h — 1. In the limit h — > the result 
for A and thus the escape rate T diverges, because the 
saddle becomes infinitely wide. This divergence is, how- 
ever, unphysical, because Eq. ( |l3| ) requires £ ee 2/wt » 1. 

For small a, the IHD formula above fails, because the 
tacit assumption that the magnetic moments approach- 
ing the barrier from the depth of the well are in thermal 
equilibrium is violated. In the LD limit, the situation be- 
comes completely different: The function g changes in a 
narrow region along the lines corresponding to the saddle- 
point energy e c , i.e., across the lines AO A' and BOB 1 in 
Fig. [I]. Diffusion in the energy space becomes very slow, 
and thus g equilibrates along the lines e = const. That 
is, C of Eq. (£>|) can be approximated as 

C(a^) = C(e), e = e(x,<p)=H(x,<p)/T. (15) 

Then the conservative part of Eq. (^) for £ vanishes. The 
rest can be averaged over the phase variable, i.e., over the 
constant-energy line, and written 
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where 
A(e) 

B(e) 




(17) 



with the averaging over conservative trajectories. 

Usually A/B ~ T/Echm, with E c hai a characteristic 
energy, can be neglected at low temperatures, and Eq. 
( |l6| ) is easily integrated. Calculating the flux across the 
part of the separatrix e = e c around the first well yields 



1 de 

1 — x 2 dip 



dx 



(18) 



i.e., A = 5/2, where S — S/T and S is the energy dis- 
sipated over the period of the precession at the saddle- 
point energy in the low-damping case, coinciding with 
the result for particles in the LD limit 11 1, the factor 
1/2 arising because of the two symmetric wells. 

If the transverse field satisfies h <C 1, then Eq. 
( |l8| ) greatly simplifies. One can set 1 — x 2 =>■ 1 ev- 
erywhere and retain only the leading term containing 
de/dx = —lax, where along the separatrix x — x(tp) = 
-2-\/hsm.((p/2). This results in 



-y/2h(l — cos<p) 



2tt 



A = ao~l x((p)dip = Saavh, 
'o 



(19) 



vanishing in the limit h — ► 0, instead of approaching Eq. 
(|l2|), and becoming invalid at £ < 1 (see below). 

Let us now study crossovers to the uniaxial regime from 
the HD and LD regimes, starting from the HD one. For 
large values of cr one expects a crossover to the uniaxial 
case for h ~ 1/cr, which is so small that the function £ 
does not yet deviate from its uniaxial form written in the 
line above Eq. ( |l2] ) . Under this assumption which will be 
checked shortly, the only effect of the transverse field is 
that in Eq. (0) / eq cx exp(£ cos c^). Integration over ip 
yields 



A = 2ira\/o-/ircxp(— £)2o(0> 



(20) 



which interpolates between Eq. ([L2|) and the third line of 
Eq. (Q), so describing a crossover between the uniaxial 
and HD regimes. If a <C 1, then at higher fields, h > a 2 , 
a crossover to the ID regime with A — 1 occurs, which is 
described by Eq. @. Thus Eq. @ applies if 1/cr < a 2 , 
i.e., a = d-^/cr ^> 1. 

Let us now consider LD for h -C 1 more accurately. 
Here A and B of Eq. |L7]) simplify to A(e) ^ (d 2 e/dx 2 ) 
-2cr and B(e) ^ ((de/dx) 2 ) = -4cre, so that Eq. |[|) 
becomes 



de 2 
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2s 



de 



dC 

de 



C 



(21) 



where the constant C is determined by the condition 



2C 



de 



exp(e - e c ) = 1, 



(22) 



Note that for £ < 1 the term l/(2e) in Eq. @ cannot 
be dropped, since — e ~ — e c = £. This is the difference 
from the standard LD case above. 
Thus, instead of Eq. ©, 



A^aC lim fl^fflfr = toayfiftfMQ, (23) 



where 

and 
/(C) 



— £ 



(24) 



dter 1 



(eWcV?) • (25) 




One has Q = 1 for f = and Q = 2 3 / 2 /tt rs 0.900 for 
£; 7^ 0, whereas the asymptotes of the function /(£) are 



(26) 



Thus for £ = the uniaxial limit, Eq. (|l2|), is recov- 
ered. In the region £ ~ 1, the function /(£) describes the 
crossover to the standard LD result of Eq. ([l9|) . The dis- 
continuous form of Q above shows that our treatment was 
not accurate enough to describe yet another crossover at 
£ ~ s/a, where a = a^fd is small in the LD case (see 
below). 

It follows that the parameters governing the uniaxial 
crossover are 



i = 2ha = HM S V/T, a = a^fd. 



(27) 



For a ^> 1, in the relevant region £ ~ 1 one has h <C 1. 
Thus in Eq. ( p^| ) the terms h cos if, adQ/dip, and d 2 (,/dip 2 
can be neglected, the derivatives with respect to y> being 
much smaller than those with respect to x. In contrast, 
the terms hsinipdC/dx and d(/dip become relevant for 
small a and should be retained. The resulting equation 
can be cast into the scaled form 



d 2 C <X 
ox ox 



2 ~<K 
dip 



d( 



(28) 







where x = \fax and the boundary conditions are C 
for x = -co and £ = 1 for x = oo . 

The "phase diagram" of the regimes for the escape 
rate of the uniaxial model with transverse field is shown 
in Fig. ||. The ID-HD crossover is described by Eq. jl^ ) 
and occurs at a 2 ~ h, i.e., a ~ VC- The ID-LD crossover 
occurs if in Eq. ( |l9| ) A ~ 1, which amounts to a ~ 1/V?- 
This crossover is described by Melnikov's formula p2[. 
The HD-Uniaxial crossover, which is described by Eq. 
(P0|), occurs for a ^> 1 at £ ~ 1. For q < 1, there 
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FIG. 2. Different regimes for the escape rate of the uni- 
axial spin system in the transverse field (a 1, h <C 1). 



are two crossovers between the LD and uniaxial regimes. 
One of them occurs at £ ~ 1 and is described by Eqs. 
( p3| ) and (|25|). Another one occurs at a ~ £ 2 . 

The latter follows from the perturbative solution for 
the escape rate at small transverse fields for arbitrary 
values of a. To second order in h, the result has the form 

T(h)/T(p) = 1 + (£ 2 /4)F(a) = 1 + h 2 a 2 F(a), (29) 



where 



F(a) = 1 



a 2 jo 



du exp 



ln(l 



2a 2 



= 1 + 2(2a 2 e) 1 /( 2Q2 ) 7 (l + l/(2a 2 ), l/(2a 2 )) , (30) 

and 7(0, z) — J Q dtt a ~ 1 e~ t is an incomplete gamma func- 
tion. The limiting forms of the above expression are 



F 



1 + 1/a- l/(2a) 2 + a»l 
\pK Jo, — 1/3 + \fna/6 + . . . , a«l. 



(31) 



The last formula shows that for a <C 1 the escape rate 
T(h) essentially deviates from its uniaxial value if a ~ £ 2 . 
This defines the other crossover mentioned above. 

The results of the numerical calculation of the thermal 
activation rate as the lowest eigenvalue of the Fokker- 
Planck equation j|] are shown in Fig. [|. For large a 
the agreement with the HD-Uniaxial crossover formula, 
Eq. (p0|), is rather good. For small a the field interval 
is restricted due to convergence problems. Nevertheless, 
there good accord with the LD crossover formula, Eq. 
(p3|), in the region £ < 1, where the standard LD result 
of Eq. JlS3| ) is invalid. 

We have shown how different damping-dependent 
regimes of thermal activation for uniaxial magnetic par- 
ticles with transverse field merge into the single uniaxial 
regime when the field tends to zero, and presented the 
complete "phase diagram" of the different regimes. The 
uniaxial characteristics appear for £ < 1, i.e., for h < l/u 
[see Eq. (||)] . It should be noted also that the transition 
from classical to quantum regimes of the escape with de- 
creasing temperature is strongly modified by proceeding 
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FIG. 3. Transverse-field dependence of the inverse relax- 
ation rate for a = 10 and different values of the damping 
constant a. Solid lines represent Eq. ([fl]) with A given by 
Eqs. © and with Q = 2 3/2 /tt » 0.900. 



to the uniaxial limit |l3[. The latter is, however, a more 
pronounced effect and it already occurs for h < 1/4 for 
our model. 
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